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Abstract. We study a time-dependent hydrogen ionization in the atmosphere of SN 1987 A during the first month after the ex- 
plosion. The model includes kinetics of hydrogen ionization and excitation, molecular hydrogen kinetics, and a time-dependent 
energy balance. The primary strong effect of the time-dependent ionization is the enhanced hydrogen ionization compared to 
the steady-state model. The time-dependent ionization provides a sufficient population of excited hydrogen levels to account 
for the observed Ha without invoking the external ^^Ni. We find that the Ba II 6142 A line in SN 1987A can be reproduced 
for the LMC barium abundance. This resolves the long-standing problem of the unacceptably high barium overabundance in 
SN 1987 A. The key missing factor that should be blamed for the "barium problem" is the time-dependent ionization. The 
modelling of the Ho- profile on day 4.64 indicates the ratio of the kinetic energy to the ejected mass » 0.83 x lO^^ergMg . 
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1. Introduction 

Spectra of type IIP supemovae (SN) combined with photo- 
metric data are the primary source of our knowledge about 
the mass, energy, chemical composition, mixing, asymmetry 
of ejecta that are of a vital importance for the verification of 
the explosion mechanism. Generally, to recover the informa- 
tion imprinted in the spectra one needs to use rather elaborated 
models. Unfortunately, the ultimate model with all the physics 
included is beyond reach, so some simplifications are unavoid- 
able. One of the assumptions accepted in the atmosphere mod- 
els is the steady-state approximation for the ionization kinetics. 

Yet, that the time-dependent effect of the hydrogen ioniza- 
tion may play a role in the SN envelope has been recognized 



by Kirshner & Kwan J1975> : they applied this effect to ac- 
count for the high Ha luminosity of the type II SN 1970G. The 
time-dependent ionization has been exploited for SN 1987 A 
to account for the high excitation of hydrogen in the outer at- 
mosphere (v > 7000 kms"') during the first 40 days (Chugai 
[1991a-). The crucial role of the time-dependent ionization at 
the late nebular epoch of SN 1987 A has been emphasized by 
Clayton et al. (1952^ and Fransson & Kozma (1993 ). 

Recently, using a time-dependent model of the hydrogen 
ionization with the kinetics of atomic and molecular hydrogen 
we have found that time-dependent effects in SN 1987 A were 
strong during the first month and that they could provide the hy- 
drogen excitation required to account for the Ha line strength 
(Utrobin & ChugailOOZ) which was underproduced in steady- 
state models of SN 1987A atmosphere. However, our previous 
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model has ignored a time-dependent energy balance; instead 
simple laws have been used to mimic the evolution of the elec- 
tron temperature in the atmosphere. 

Keeping in mind the importance of the non-steady state ion- 
ization effects for the understanding of the phenomena seen in 
early SN 1987 A, we revisit this issue. To this end we upgrade 
the atmosphere model that apart from the time-dependent ki- 
netics includes also the time-dependent energy balance. Our 
main goal is to compute the ionization and excitation of hydro- 
gen in the atmosphere of SN 1987 A for different epochs, and to 
compare the calculated Ha line with the observations. We also 
address a highly intriguing problem of tremendous Ba over- 
abundance that has been recovered by several studies (Mazzali 
et al. 1992 and references there). Our preliminary guess was 
that the Ba overabundance problem was caused by the strong 
underestimation of the Ba II ionization fraction. 

We begin with a question, why time-dependent effects are 
important (Sect. |2} and then describe in detail the proposed 
model (Sect. The Ha strength on day 4.64 is used to con- 
strain the kinetic energy of the hydrodynamic model Sect. 14.1! 
the computations of the ionization and temperature structure of 
the atmosphere along with the Ha profiles for selected phases 
are presented in Sect. 14.2! while in Sect. 14.31 we address the 
problem of high Ba overabundance. 

We adopt here the radial velocity of SN 1987 A of 286.5 
km s ' (Meaburn et al. il995i and the distance of 50 kpc. 
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2. Why time-dependent? 

A widely shared view that the steady-state statistical equilib- 
rium is a reasonable approximation for SN IIP atmospheres re- 
lies on the claim that the recombination time is much smaller 
than the SN age 

(aN^r^ « t , (1) 

where a is the coefficient of recombination to the excited states, 
A^e is the electron number density. Adopting a ~ 3 x 10"'^ 
cm^ s"' and a typical SN age at the photospheric epoch f ~ 
10^ - 10^ s one finds that the requirement Q is satisfied 
for A^e » 10^ cm"^. This condition is actually met at the 
photosphere thus supporting the steady-state approximation. 
Meanwhile, it is violated in the outer, high velocity layers, 
where Ng may drop as low as 10^ cm""*. Therefore, generally, 
the time dependent approach is needed at least, if one addresses 
the situation in the outer layers. 

Moreover, it turns out that the steady-state approximation is 
violated also in the inner atmosphere, where the condition Q 
holds. The point is that this inequality has little in common with 
the condition of the steady state for multilevel atom irradiated 
by the optical continuum, because of the efficient reionization 
of recombined atom. The importance of this phenomenon has 
been first recognized for the problem of the hydrogen recom- 
bination in the early universe, where La trapping along with 
the reionization from the second level results in the signifi- 
cant increase of the recombination time fi-ec - -Ne/(dNs/dt) 
(Zeldovich etal. lT968t . 

The net recombination rate dNg/dt refers to transitions to 
the ground state that either do not emit La or followed by the 
escape of the emitted La quantum without resonance scatter- 
ing. Let the probability of such a transition is W2i- Then in the 
"two level plus continuum" approximation the net recombina- 
tion rate (neglecting the expansion effect) is 

dN^ , 

— ^ = -aN^,W2i . (2) 
dt 

With all the major processes from the second level taken into 
account the probability W2\ is 

A2q + A^e?21 +A2ly821 

^21 = -. 77 -. S~ ' (3) 

P2 + A2q + Aeg'21 + A21/321 

where P2 is the photoionization rate from the second level; A2q 
is the probability of the two-photon decay evaluated for the 
equipartition of 2s and 2p levels; q2i is the collisional deex- 
citation coefficient; A21J621 is the La escape probability. 

Let us estimate W21 for the neutral hydrogen number den- 
sity Ni > 10^ cm"-' and the electron number density Ae < 10** 
cm when the two-photon transition dominates over colli- 
sional transition and La escape, i.e., W2\ ~ A2q/P2- This sit- 
uation is typical for the inner atmosphere of SN IIP. With the 
photospheric temperature of T x 5000 K one gets P2 ~ 10*W 
s where W is the dilution factor. Close to the photosphere 
(for W > 0.1) one finds W21 < 2 x lO""*. The effective recombi- 
nation time is then 

?rec = -NJ(dNJdt) = l/(aAeW2i) ~ 500/(a'Ae) . (4) 




Ig A (A) 



Fig. 1. The combined UV (lUE) (Pun et al. 1995^ and opti- 
cal (CTIO) (Phillips et al. .1988J spectra of SN 1987A (thin 
solid line), the black-body continuum at the calculated effec- 
tive temperature of the high-energy model (dotted line), and 
the observed continuum (thick solid line) between days 4.64 
and 44.2 L 

This shows that the effective recombination time is by a factor 
of l/vv2i ~ 500 greater than the recombination time to excited 
states l/(Q'Ae). For Ae ~ 10*^ cm"-' one obtains free > 10^ s 
> f. We thus conclude that the steady-state approximation may 
break down in the dense inner atmosphere as well. 

These qualitative considerations demonstrate that the hy- 
drogen recombination in the atmosphere of SN IIP at the pho- 
tospheric epoch is essentially a time-dependent phenomenon. 

3. The model and input physics 

A self-consistent treatment of the SN atmosphere requires the 
hydrodynamic modelling with time-dependent radiation trans- 
port, energy balance, and ionization kinetics. At present such 
a general approach is beyond reach of computation possibil- 
ities. To study time-dependent effects in the atmosphere of 
SN 1987 A and yet to make the problem solvable, we sacrifice 
the radiation transfer in the atmosphere using, instead, a simple 
description of the continuum radiation field. 

3.1. Continuum radiation 

The SN atmosphere is very opaque in the Lyman continuum 
(v > vi), so the photoionization by the photospheric Lyman 
continuum can be neglected. On the other hand, the hydrogen 
ionization by the recombination Lyman continuum in the at- 
mosphere can be easily taken into account via cancelling the 
photoionization of hydrogen from the ground level and the re- 
combination to the ground level in the rate equations. In the at- 
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mosphere the radiation field in continuum at these frequencies 
is considered as being in the thermal equilibrium with matter 

The important simplification we adopt is that the continuum 
in the range of v < vj corresponds to the free-streaming case, 
i.e., the average intensity of the continuum in the atmosphere is 
Jv - Wly, where /y is the specific intensity of the photospheric 
radiation. The photospheric radius R^^ and effective tempera- 
ture Tett are provided by the hydrodynamic model. 

The free-streaming approximation is a good approximation 
in the visual band, where the optical depth of the atmosphere is 
quite low. However, for the ultraviolet (UV) radiation between 
Lyman and Balmer edges this approximation is indeed rather 
artificial, because of the large optical depth due to the UV line 
opacity. As a result, the average intensity of the UV radiation in 
the atmosphere generally should differ from the free-streaming 
approximation. 

Nevertheless, this approximation is not as bad as one would 
think at first glance. It predicts behavior of 7 oc r"^ in the 
most of the atmosphere. In the case of pure scattering with 
the constant scattering coefficient ixir) - const) one obtains 
for the extended atmosphere 7 oc [1 -h \.5Rx(Rlr - 1)] where 
R is the outer radius of the atmosphere (Chandrasekhar'1934 l. 
According to this expression the behavior of J{r) is steeper than 
r"^ for outer layers (r ~ R) and J{r) oc r ' for r ^ R. So, our 
approximation corresponds to some intermediate range in the 
atmosphere. Moreover, since the scattering of UV metal lines 
is liable to photon splitting, the interaction of the UV radiation 
with matter will be close to the absorption rather than scatter- 
ing. We expect, therefore, more steep behavior of J{r) in the 
inner layer compared to the pure scattering case. In that case 
the law J oc r^^ may be a sensible approximation in the deeper 
layers as well. 

To further allow for the uncertainty of our approximation 
we consider two extreme cases: (i) the photospheric brightness 
is black-body with the effective temperature (model A); (ii) the 
photospheric brightness corresponds to the observed spectrum 
(model B) (Fig. 0. The real situation is somewhere between 
these two cases. 



3.2. Hydrodynamic model 

The density distribution, chemical composition, radius of the 
photosphere and effective temperature are provided by the hy- 
drodynamic model of SN 1987A (see for details Utrobin 2004 1. 
We consider the non-evolutionary model N of Utrobin ( 2004 ) 
with 18 Mq ejecta, kinetic energy of 10^' erg, and ^^Ni mass of 
0.073 Mq. For the distance modulus to the LMC of m - M - 
18.5 mag, color excess of E{B - V) = 0.15 mag, and an in- 
terstellar extinction of Ay - 3.1E{B - V) mag a progenitor of 
SN 1987A, Sanduleak -69°202, had the radius of 46.8 Rq. 
To produce the presupernova density distribution the evolu- 
tionary model 120n2ae of Woosley et al. il997> was rescaled 
to adopted radius and ejecta mass. We assumed for the outer 
layers the circumstellar helium abundance Y - 0.441 (Wang 
119911 Lundqvist & Fransson I1996> and the LMC metallicity 
Z - 0.004 (Dufour"1984'l. The structure of the presupernova 
outer layers was modified to increase density compared with 
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Fig. 2. Comparison of the low and high-energy models of 
SN 1987 A. Top: (a) bolometric light curve of the low-energy 
model (thick solid line) compared with the observations of 
SN 1987A obtained by Catchpole et al. ( 1987"1988"l (open cir- 
cles) and Hamuy et al. (.1988 ) (open triangles), and (b) the Ha 
Une profile on day 4.64 computed on the basis of this hydrody- 
namic model for the black-body continuum (model A) (dotted 
line) and the observed continuum (model B) (thick solid line). 
The calculated spectra are overplotted on the observed spec- 
trum (Phillips et al. '1988'^ (thin solid line). Bottom: (c) bolo- 
metric light curve of the high-energy model and (d) the corre- 
sponding Ha profile on day 4.64. 



the evolutionary model of non-rotating single star This cor- 
rection was needed to fit the observed bolometric light curve 
of SN 1987A (Fig. 2^) with a moderate mixing of radioactive 
^^Ni in the velocity range < 2500 kms ' (Utrobin EUO^t . The 
final model is dubbed as the low-energy model. 

We consider also the high-energy version of the above 
model. This model with the kinetic energy 1.5 x 10^' erg dif- 
fers from the low-energy model by the presupernova radius of 
35.0 Rq, ^^Ni mass of 0.075 Mq, and the adopted color excess 
of E(B - y) = 0.17 mag (Michael et al. 2003 1. The chemical 
composition of the outer layers is typical to the LMC chemical 
composition, X = 0.743, Y = 0.251, and Z = 0.006 (Dufour 
1984). This model agrees with the observed bolometric light 
curve of SN 1987 A for ^^Ni mixed in the velocity range < 3000 
kms-i (Fig.El;). 

3.3. Line transfer 

Starting the first day the envelope of SN 1987 A expands ho- 
mologously (Utrobin .2004 so the Sobolev approximation is 
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applicable. In this approximation the line transfer is described 
by the escape probability (Sobolev [T960l Castor ll970t 



Pul 

where 



1 - exp{-Tiu) 



1X6 I Si 

Ttu = flu^lut\Nl-^N, 



(5) 



(6) 



is the Sobolev optical depth in a transition from level / to 
level u. Here Ni, gi, Nu, g„ are number densities and statistical 
weights of atoms in the lower and upper levels, respectively, //„ 
is the oscillator strength, and /l;„ is the line wavelength. In the 
Sobolev approximation the frequency-averaged mean intensity 
of the line is 



where 



giNu 



(7) 



(8) 



is the line source function and v„/ is the line frequency. The net 
frequency-integrated line emissivity is 



J]ui = —hviuNuAuiPul 1 — 

4n \ Siu 



(9) 



where A„/ is the Einstein spontaneous emission probability for 
the transition from level u to level I. 

For the very opaque La line the Sobolev approximation is 
violated, if one consider the ha radiation transfer, but remains 
valid, if one is interested in the computations of the escape rate 
(ChugaidH. 

3.4. Gamma-rays and positrons deposition 

The gamma rays with energy of about 1 MeV from the de- 
cay chain ^^Ni ^^^Co — >^''Fe deposit their energy through 
Compton scattering by free and bound electrons. The Compton 
electrons lose their energy through the Coulomb heating of free 
electrons and the ionization and excitation of atoms and ions. 
The rates of the heating, non-thermal ionization and excitation 
for atoms and ions are taken according to Kozma & Fransson 



JI992> . The gamma-rays transport is calculated in the approx- 
imation of the effective absorption opacity of 0.06 Fe cm^ g"', 
where is the number of electrons per baryon. Positrons are 
assumed to deposit their energy locally. 

3.5. Rate equations 

The following elements and molecules are calculated in time- 
dependent non-LTE chemical kinetics: H, He, C, N, O, Ne, Na, 
Mg, Si, S, Ar, Ca, Fe, Ba, H , H2, H+, and H+. All elements but 
H are treated with the three ionization stages. Neutral hydrogen 
and ions Mg II, Fe II, and Ba II are modelled with 15, 9, 30, and 
17 levels, respectively, while other atoms and ions are assumed 
to consist of the ground state and continuum. The reaction net- 
work includes all bound-bound and bound-free, radiative and 



collisional processes for the atoms and ions with a detailed 
level structure (Appendices lA. lHA.4t . all bound-free, radia- 
tive and collisional processes for the two-state atoms and ions 
(Appendix IA.5> . We consider kinetics of hydrogen molecules 
that may affect the hydrogen ionization degree. Seven radiative 
and 37 collisional processes for the molecules are taken into 
account (Appendices lB. ltfB.2> . 

For one-dimensional, spherically symmetric flow the equa- 
tion of continuity simplifies to 



dp 

dt 



1 



2 + 



dlnv\v 
dlnr 



(10) 



^exp ^exp 

where d/dt is the Lagrangean time derivative; p is the density 
of the matter; r and v are radius and velocity, respectively. The 
net rate of transitions between level k and all other levels of 
neutral atom z by atomic processes is 



dt 



Jill 

Sku 



- N-. 



1 - 



- N^».kPk 



u>k 

l<k 

+ A^e ^ {N^o,uqi,k - N^o^i.qku) + N^^N^qck 

u>k 

- A'e ^ [NjO^kqu - N^ajqikj - NeN^o^kqk,c 



Jill 

Sik 



l<k 



k,c 



(11) 



where N^o k is the number density of neutral atom z at the level 
k; Nz+ is the number density of singly ionized element z; A^e 
is the number density of free electrons; 7/„ is the frequency- 
averaged mean intensity of the line 0; 5/,, is the line source 
function (|S}; ak and Pk are the total radiative recombination 
coefficient, including the spontaneous and induced processes, 
and the photoionization rate of neutral atom z in level k, re- 
spectively (Mihalas 1978 1; qiu is the electron collisional ex- 
citation rate from level I to level u; qui - {NiO ilN^a ^yqu, is 
the corresponding deexcitation rate and the ratio (N^aj/N-iP i,)* 
for the local thermodynamic equilibrium (LTE) is given by the 
Boltzmann formula; qk^c is the electron collisional ionization 
rate from level k; q^^k - (N^o^k/Nz+Neyqk^c is the correspond- 
ing recombination coefficient and the LTE ratio {NzP^k/Nz+Ne)* 
is given by the Saha equation; and R^^ are the non-thermal 
excitation and ionization rates of neutral atom z in level k, re- 
spectively (Kozma & Fransson 1992 1. 

The rate equation for the number density A^h»,2 of neutral 
hydrogen in level 2 is determined by the net rate of flow of par- 
ticles into the unit volume with a characteristic time t^xp (I10> . 
the net rate (II 1> . hydrogen two-photon decays, and other de- 
population processes: 



H»,2 



dt 



^exp 



dNnO ; 



dt 



- ksNwNw^is - keNv«Nii«^2p 

- ksNiiO iNuo,2s - kgNii0jNii0,2p 

- kuNii^Nyp^2 , 



(12) 



where A^h",2s - 0.25A^h°,2 ™d A^H°,2p = 0.75A^h°,2 the num- 
ber density of hydrogen levels 2s and 2p calculated assuming 
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only collisional transition between these States; A^H" = Z A^h"* 

k=l 

is the number density of neutral hydrogen; A^h, is the number 
density of H2; Aj^j^ is the probability of the hydrogen two- 
photon decay; k^, k^, k^, kg, and ki4 are the rate coefficients 
for the corresponding reactions from Table IB. l1 The rate equa- 
tions for the number density of neutral hydrogen in level 3 and 
higher levels k > 3 are 



H°,3 



H°,3 



dt 

and 

dt 



dN, 



H»,3 



^exp 

^uA^HjA^ho.s 



dt 



- ^vA^HiA^H" ; 



A^HOJ 



^exp 



dNiioj 



dt 



■ kuNH2NH<>j 



(13) 



(14) 



respectively. Similarly, the rate equation for the number density 
Nz+,k of ions = Mg^, Fe^, Ba^ in level k simply is 



dN^ 
dt 



^exp 



dN,., 



dt 



(15) 



The rate equations for ionized hydrogen, negative hydrogen 
ion, and molecular hydrogen H2, H2 and H3 are given in 
Appendixld 

For the two-state atoms and ions the net rate of transitions, 
for example, from neutral atom z into ion z^ by ionization and 
recombination processes is given by 



dN,* 



dt 



NAPz^+N^qzo.c + RzoJ 



(16) 



Using the net rate ( I16t the rate equations for elements with 
three ionization stages are easily generalized: 



(17) 



dN,. 


-_^ + 

texp 


dN,. 




dN,2. 


dt 


dt 




dt 



and 

dN,2. 



dt 



dN,2. 



dt 



(18) 



In particular, the rate equation for singly ionized oxygen taking 
account of the charge transfer processes in collisions of with 
O and 0+ with H is 



dNo* 


- _^o* + 


dNo*' 




'dNo2.' 


dt 


texp 


dt 




dt 



+ knoNn^Noo - kouNo^Nif, , 



(19) 



where A^h+ is the number density of ionized hydrogen. The net 
rate J16> for the ions with a detailed level structure (Mg^, Fe^, 
and Ba^) is calculated including all relevant atomic processes. 

To close the system of the above rate equations J12H15I 
IC.1HC.5I [T7HT9t we use the particle conservation for hydro- 
gen 



Nuo + Nii^ + Nii- + 2Nu, + 2NHt + 3Nn* = 



(20) 



and for other elements from He to Ba with three ionization 
stages 



muAz 

and the charge conservation 



A^H^ - A^H- + A^m +Ny^ + Y^ (N,. + 2N,2.) = A^e 



(21) 



(22) 



z=He 



where is the mass fraction of element z, and Az is its atomic 
weight. 

3.6. Gas energy equation 

The equation of state of a perfect gas can be written as 



(N, + N,)kT, 



(23) 



where Pg is gas pressure and Nt is the total number density of 
atoms, molecules, and ions of all types. The specific internal 
energy of particles (per gram) is the sum of the translational, 
excitation, and ionization energies: 



p - __! IkT + F + F- 

2 p 

where 

Ba+ 



Eexc ~ 

z=H k>\ 

is the excitation energy and 

Eion = /hA^h+ - /h-A^h- + Ih.Nhi + /h3A^h 

Ba 
z=He 



(24) 



(25) 



(26) 



is the ionization energy. Here e2,k is the excitation energy of 
state k for the corresponding atoms and ions, and /h, Ih-, Ih^, 
/h,, Iz", h+ are the ionization potentials of the corresponding 
atoms, molecules, and ions. 

The gas energy equation, the first law of thermodynamics 
for the material, including radiative losses in continuum and 
lines, the Compton cooling, and non-thermal heating thus reads 

dE, ^ p 5 /1\ An , 



An 
P 



/'Ba+ \ 
>,z=H l,u 



+ E' 



where 



c AkT^ I T, 



(27) 



(28) 



is the net emissivity of the Compton cooling (Weymann'19661. 
Here Ky and ril are the monochromatic true absorption coef- 
ficient corrected for stimulated emission and the monochro- 
matic thermal emissivity, including bound-free and free-free 
processes for all atoms and ions, respectively (Mihalas ll978> : 
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rjy is the two-photon emissivity of hydrogen; ?/„/ is the net 
frequency-averaged emissivity of line (|9jl; is the rate of en- 
ergy deposition from the radioactive decays; cr^ is the Thomson 
scattering cross section; J is the integrated mean intensity of 
continuum; - (nJ/crj^y^'^ is the radiation temperature and 
ctr is the Stefan-Boltzmann constant. The energy equation de- 
termines the electron temperature in the SN atmosphere which 
is the same for neutrals and ions. 



3.7. Computational method 

A hydrodynamic model with 300 mass zones provides us the 
velocity v(r, f) and density p(r, t) profiles, the photospheric ra- 
dius /?ph(f) and effective temperature Tef[(t). The simplified 
description of the continuum radiation field as diluted photo- 
spheric radiation permits us to solve for a given mass zone the 
rate equations ( fnUTSl IC.1HC.5I [T7HT9t . the particle conserva- 
tion equations ( 1201 12 1> . the charge conservation equation i22\ . 
and the gas energy equation (|^} independently of other mass 
zones. 

As the supernova envelope expands the photosphere moves 
inward the envelope. In a steady state the physical values in 
the atmosphere at certain moment t are determined by the in- 
stant density and radiation field. In a time-dependent approach 
one has to start from a moment fph at which a given mass zone 
crosses the photosphere, and then to follow the time evolution 
until t. To obtain the distribution of physical values in the at- 
mosphere this procedure should be applied to each mass zone. 

The simultaneous solution of the above equations (I12H14I 
IC.1HC.5| I17H22I I27> is achieved by reducing them to a system 
of ordinary diff'erential equations. In particular, the gas energy 
equation ( I27> is reduced to the ordinary differential equation 
for electron temperature Tg. To properly formulate the prob- 
lem of time evolution of physical values for a given mass zone, 
initial data must be specified at the photospheric level. This is 
done by solving the equations ((T2HI3I IC.1HC.5I fTTti^ at the 
photospheric level assuming the steady state for a given elec- 
tron temperature Tg. This is the simplest way to model the non- 
equilibrium conditions at the SN photosphere. The derived sys- 
tem of ordinary differential equations is stiff and is integrated 
by the implicit method of Gear (1971 1 with an automatic choice 
of both the time integration step and the order of accuracy of 
the method. 

The recovered time-dependent structure of the atmosphere 
is used then to calculate spectra at selected epochs. The spec- 
tra are modelled by Monte Carlo technique with relativistic ef- 
fects included (Mihalas lTWSl Jefferv ll993H1995> . We suggest 
that the photosphere diffusively reflects the incident photons 
(Chugai & Utrobin 2000 l. The line scattering is generally non- 
conservative and is described in terms of the line scattering 
albedo. The Thomson scattering on free electrons and Rayleigh 
scattering on neutral hydrogen are taken into account. 



4. Results 

4.1. Early Ha and the hydrodynamic model 

Preliminary modelling of Ha has revealed that at the early 
epoch the Ha line strength was sensitive to the density of the 
outer layers. A similar sensitivity of the Ho- strength to the den- 
sity was found in the steady-state model (Eastman & Kirshner 
.1989) . We use therefore the Ha on day 4.64 to constrain the 
kinetic energy of the hydrodynamic model. 

Assuming the initial electron temperature equal to the ef- 
fective temperature Tg = Teff we solved the system of rate 
and energy equations and then computed Ha profile for the 
low-energy model, i.e., for the model N with the kinetic en- 
ergy E - 10^' erg (Utrobin 2004). The corresponding bolo- 
metric light curve and the calculated Ha profile on day 4.64 
are shown in Figs.|5^ andl^'. The computed Ha line has too 
weak absorption at high velocities compared to the observa- 
tions (Phillips et al. ll988> . The situation, however, is markedly 
improved for the high-energy model with £ = 1.5 x 10^' erg 
(Figs.|2j; and|2ll). The larger strength of the Ha absorption in 
the high-energy model is the direct outcome of the higher den- 
sity in the outer layers. Thus, the strength of the Ha absorp- 
tion at high velocities indicates the favoring energy-to-mass ra- 
tio E/M X 0.83 X 10^° erg Mg'. We adopt the model with the 
high energy as the standard model for our computations of the 
SN 1987A atmosphere. 

4.2. Evolution of hydrogen ionization and Ha 

With the adopted photospheric continuum and the initial elec- 
tron temperature we are able to compute for given moments the 
distribution of the essential physical values in the atmosphere, 
including the hydrogen level populations. The initial electron 
temperature at the photosphere is adopted equal to the effec- 
tive temperature T^ff which is taken from the hydrodynamical 
calculations (Utrobin . 2004j . 

4.2.1. Day 4.64 

On day 4.64 the Ha profile computed in the time-dependent 
model A fits the observed one fairly well (Fig. |3^), while the 
model B underproduces the emission component (Fig. |3j>)- 
This suggests that the average intensity of the UV radiation in 
the atmosphere at this and earlier epochs is closer to that of the 
model A. The steady-state model A produces rather strong Ha 
line although both the absorption and emission components are 
markedly weaker than in the observed line, while in the steady- 
state model B the Ha line is extremely weak. The weakness of 
the Ha line in the steady-state model reflects the simple fact 
that the steady-state ionization is lower than that in the time- 
dependent model. 

As expected, the electron temperature in the model A is 
larger than in the model B since the radiative heating rate in the 
model A with the black-body UV continuum is higher than in 
the model B (Figs.|3l; andlsjl)- The remarkable property of the 
distribution of electron temperature in the atmosphere for both 
models is a fast drop from the initial value at the photosphere of 
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Fig. 3. Ha line profile and the structure of the atmosphere 
on day 4.64. Panel (a); Ha profile computed with the time- 
dependent model A (thick solid line) and with the steady-state 
model A [dotted line), and overplotted on the observed profile 
obtained by Phillips et al. C1988,) (thin solid line). Panel (b): 
the same as panel (a) but for the model B. Panel (c): the cal- 
culated electron number density (thick solid line) and electron 
temperature (thin solid line) in the time dependent model A. 
The dashed line shows the initial electron temperature Tf^ at 
the photosphere as a function of velocity, while the dotted line 
is the distribution of the electron temperature T^'^ evaluated 
with the purely adiabatic cooling. Arrow indicates the photo- 
spheric level. Panel (d): the same as panel (c) but for the model 
B. 

Teff to ~ (0.75-0.80) Teff. This drop results from a simple fact 
that the electron temperature calculated in the diluted radiation 
field is lower than the effective temperature. The same physics 
leads to a well known relation for the boundary temperature of 
gray atmosphere To 0.8iretf (see, e.g., Mihalas"1 978> . 

Although the electron temperature is relatively high in the 
most of the atmosphere, Te > 4000 K, (Figs. ^ and |3ll) it 
is insufficient to maintain the large optical depth of Ha demon- 
strated by the time-dependent models. The primary mechanism 
of the population of the second level is the recombination to ex- 
cited levels. 

Noteworthy, the populations of excited levels are small 
numbers that are controlled by the relatively fast transition 
rates. As a result the relaxation times for the level populations 
of hydrogen are small: A^h°,*/(<^-^h°,*:/'^0 ^ t- So the popula- 
tions of excited levels may be considered in the steady-state 
approximation. We conclude that the population of the sec- 
ond hydrogen level, and the Ha optical depth as well, is deter- 
mined mainly by the current values of the electron concentra- 
tion and the photoionization rates, which are essentially time- 
dependent. 



Fig. 4. Ha line profile and the structure of the atmosphere on 
day 8.69. See Fig.|3llegend for details. 



4.2.2. Day 8.69 

On day 8.69 the both A and B models reproduce satisfacto- 
rily the absorption component of Ha (Figs. |^ and|3)). The 
strength of observed emission lies between two models which 
implies that the actual radiation field in the UV band is interme- 
diate between A and B cases. The steady-state approach in both 
models predicts very weak Ha line at this epoch. This is qual- 
itatively consistent with the steady-state modelling of Mitchell 
etal. (2002.). 

The electron concentration shows a kink at about 14000 
km s"' (Figs. 13; and|3l) that was only barely seen at the previ- 
ous date (Figs.|3j; andlsjl)- The position of this kink correlates 
with the transition from the full ionization to partial ionization 
at the photosphere that occurred after f ^ 2 days. 

4.2.3. Day 14.68 and 19.68 

The calculated Ha profiles on day 14.68 (Figs.|5t and|5j)) and 
19.68 (Figs.|6^ andl^Ji) show similar behavior. The model B fits 
to the observed profile fairly well, whereas the emission com- 
ponent in the model A is unacceptably strong. This indicates 
that around these epochs the UV continuum radiation is closer 
to that of the model B. The steady-state model at this phase is 
very weak with practically zero absorption and emission com- 
ponents. 

At these epochs the adiabatic cooling in the model B dom- 
inates at the low velocities as indicated by the close run of 
temperatures Te and T^'^ (Figs. |5jj and (SJl), and the most of 
the atmosphere is cooled below ~ 2000 K. In the model A the 
heating due to the photoionization is strong enough to maintain 
rather high temperature ~ 4000 K throughout the atmosphere. 
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Fig. 5. Ha line profile and the structure of the atmosphere on 
day 14.68. Observed profile is given by Hanuschik & Dachs 
([MH). See Fig.Ellegend for details. 



Fig. 7. Ha line profile and the structure of the atmosphere on 
day 29.68. Observed profile is given by Hanuschik & Dachs 
a988 1. See Fig.|3legend for details. 
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Fig. 6. Ha line profile and the structure of the atmosphere on 
day 19.68. Observed profile is given by Hanuschik & Dachs 
(fi^SSt . See Fig.Ellegend for details. 



4.2.4. Day 29.68 

On day 29.68 the model A leads to the unacceptably strong 
Ho- emission (Fig.Et)^ while the model B reproduces the over- 
all strength of the absorption and emission components of Ho- 
(Fig-Eb)- The Ha line in the steady-state models is extremely 
weak, although slightly stronger than on day 14.68 and 19.68. 




10 

V (10^ km s-') 

Fig. 8. The Ha optical depth on day 29.68 in the model A. The 
local minimum at 5000 kms"' is a characteristic feature re- 
quired for the explanation of the blue peak. 



The behavior of electron temperature T^. and electron con- 
centration A^e is very different for both models (Figs. 0; and 
01) that emphasizes the sensitivity of the result to the behavior 
of the UV radiation field in the atmosphere. Yet, it is highly 
remarkable that even the time-dependent model B, in which 
the UV continuum intensity is possibly underestimated, repro- 
duces the strength of the Ha line in the whole range of the 
radial velocities. This means that the admixed ^^Ni in the outer 
layers is not needed to account for the Ha line in SN 1987 A at 
the age of f < 30 days. 

Between day 20 and 29 the additional blue and red peaks 
emerge in Ha profile, a phenomenon known as "Bochum 
event" (Hanuschik & Dachs 1 988 >■ We believe that the red 
peak is related with the high velocity ^^Ni clump in the far 
hemisphere, whereas the blue peak reflects a minimum in 
the spherically-symmetric radial distribution of the Ha optical 
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Fig. 9. Evolution of Ba II 6142 A line from day 8.69 till 19.68 
in SN 1987A. The Ba II 6142 A profiles, computed with the 
time-dependent models A (dotted line) and B (thick solid line), 
are overplotted on those observed by Phillips et al. ("1988 1 and 
Hanuschik & Dachs (1988. 1 (thin solid line). Left panel (a): 
models for the relative Ba overabundance of 1.5. Right panel 
(b): models with the suppressed far UV continuum. 



depth at « 5000 kms ' (Chugai [TWT51 Utrobin et al.\MB- 
The mechanism responsible for such a minimum is still not 
identified confidently. Both A and B models do not show any 
signature of the blue peak in Ha profile. Remarkably, however, 
the model A reveals an essentially non-monotonic behavior of 
the Ha optical depth T(Ha) (Fig. |S} with the position of the 
minimum at 5000 kms"', just at the place that corresponds 
to the blue peak. Unfortunately, the optical depth at this min- 
imum is too high to produce the blue peak in the Ha profile. 
The required value should be 30 times lower However the very 
possibility that the non-monotonic behavior of t(v) may arise 
in a model with monotonic evolution of the photospheric pa- 
rameters and monotonic density distribution is very promising 
for the explanation of the blue peak in the Ha profile at the 
Bochum event phase. 

Previously, using the time-dependent model we found that 
the enhanced hydrogen neutralization at about 4000 kms ' 
governed by H" and Hj might be responsible for the local min- 
imum of A^e and thus, for the minimum of the Hff optical depth 
T(Ha) (Utrobin & Chugai 2002 1. However, in that model we 
adopted rather simple behavior of electron temperature T^. in 
the atmosphere and at the photosphere. As a result the role of 
the molecular processes has been significantly exaggerated. 



4.3. Barium problem 

The "Ba problem" is the disparity between the high Ba over- 
abundance factor ^ =(Ba/Fe)g7A/(Ba/Fe)LMC ~ 5 derived from 
Ba II 6142 A line (Mazzah et al. ll992> and the theoretical pre- 
dictions of i-process calculations. The ^-process nucleosynthe- 
sis in the massive stars evolving to the presupernova star of 



SN 1987 A (Prantzos et al. ll990t yields the Ba overabundance 
of ^ ~ 1.3-1.4 in the hydrogen-rich envelope, assuming that a 
total ejecta mass is 15-18 Mq, a mass-cut is at 2 Mq (Woosley 
& Weaver 1995 1, and the synthesized Ba is completely mixed 
over the ejecta. In fact, it is highly difficult to imagine that Ba 
produced at the He-burning stage can be significantly mixed 
with the hydrogen envelope. More realistic would be a situa- 
tion with the LMC abundance of Ba in the hydrogen envelope. 

In our time-dependent approach we recover the distribution 
of the electron concentration that fits to the Ho- profile on day 
14.68 and 19.68 in the model B. Therefore we hope that our 
model can recover a more realistic abundance of Ba from Ba 
II 6142 A Une. We compute the Ba ionization and excitation, 
using the multilevel Ba II model described in Appendix IA.4I 
for the two choices of the far UV (10.0 < hv < 13.6 eV) ra- 
diation in the SN envelope. In the first version the continuum 
in this energy range is the diluted photospheric radiation — our 
standard assumption for models A and B (Fig.[0. In the second 



version we follow the suggestion of Mazzali et al. J1992> that 
the far UV continuum should be essentially suppressed by the 
line-blocking (in fact Rayleigh scattering also contributed to 
the blocking effect). This assumption essentially increases the 
ionization fraction of Ba II and thus decreases the overabun- 
dance factor We, however, modify this assumption making it 
more realistic: we adopt that the far UV (10.0 < hv < 13.6 
eV) radiation is in equilibrium with matter and is character- 
ized by the local electron temperature of our time-dependent 
model. This is the same approximation as accepted above for 
hv > 13.6 eV. Unfortunately, the lUE observations cannot pro- 
vide a confident verification of our assumption about the far 
UV flux in SN 1987 A because of the dominant contribution of 
the light from the stars 2 and 3 to this band (Pun et al. '1995V 
We find that the accepted suppression of the far UV continuum 
does not affect the computed Ha hne for A and B cases at all. 

The Ba II 6142 A line profile calculated for three phases 
preceding the Bochum event are shown in Fig. |5] In the first 
version of the continuum the Ba relative abundance of ^ = 1 .5 
is required to reproduce the Ba II 6142 A line (Fig.|9^), while 
in the second case the LMC abundance of Ba (^ = 1) is quite 
enough (Fig.|9j)). We find that for ^ = 5 the blue wing of the 
line can be reproduced in the model B on day 14.68 and 19.68 
with the photospheric far UV radiation. Note that Mazzali et al. 



( I1992> obtained a similar overabundance in the model with the 
essentially suppressed far UV radiation. In our model with the 
suppressed far UV radiation the LMC abundance of Ba suffices 
to account for the Ba II 6142 A line on day 14.68 and 19.68 
(Fig.|5J)). Note that in this case the model B produces too strong 
Ba II 6142 A line on day 8.69, while the model A produces too 
weak line. This is consistent with our previous conclusion that 
at this epoch the Ha line requires the intermediate continuum 
between model A and B. 
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We come to the remarkable result: the model B with the 
suppressed far UV radiation that was also assumed in the model 
of Mazzali et al. (1992i recovers the LMC abundance of Ba 
from the Ba II 6142 A line. Therefore, we conclude that at 
present there is no serious reason anymore to suggest the over- 
abundance of Ba in the hydrogen-rich envelope of SN 1987A. 
The "Ba problem" in SN 1987 A was mainly a by-product of 
the steady-state approximation used for the determination of 
hydrogen ionization. 



5. Discussion and conclusions 

Our goal has been to explore the role of time-dependent effects 
in the ionization and excitation of hydrogen in the atmosphere 
of SN 1987 A. We have developed the model that calculated the 
time-dependent kinetics of relevant processes together with the 
time-dependent energy balance. The computations have con- 
firmed our previous conclusion (Utrobin & Chugai I2002> that 
the effect of the time-dependent ionization in SN 1987A at the 
photospheric epoch is strong: the hydrogen ionization in the at- 
mosphere is significantly larger and, as a result, the hydrogen 
excitation is higher as well. The Ha formation at the photo- 
spheric epoch during the first month is essentially controlled 
by the time-dependent ionization effect and the additional ex- 
citation by the external ^^Ni is not needed to account for the 
strength of Ha line. 

Although we have failed to reproduce the additional blue 
peak in the Ha line (Bochum event) that emerged between day 
20 and 29 (Hanuschik & Dachs 1988 1, we have found that on 
day 29.68 one of the models revealed a non-monotonic behav- 
ior of the Ha optical depth with the positions of maximum and 
minimum that were suggested by the Ha profile at this epoch. 
This result raises a hope that the blue peak will be explained 
in a more advanced time-dependent model with a correct treat- 
ment of the UV radiation transfer 

At the epoch preceding the Bochum event the time- 
dependent model produces the realistic distribution of the 
hydrogen fractional ionization in the atmosphere. With this 
structure we are able to reproduce the Ba 11 6142 A line in 
SN 1987A on day 14.68 and 19.68 with the LMC abundance of 
Ba. We believe, therefore, that the Ba problem, i.e., unaccept- 
ably high Ba overabundance in SN 1987 A, is resolved. The 
missing physics of the former models is identified: the time- 
dependent ionization that provides a higher electron concentra- 
tion in the atmosphere and thus more efficient recombination 
of Ba III into Ba II compared to the steady-state models. 
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Appendix A: Atomic data 

A.1. Hydrogen 

For hydrogen we use a model atom consisting of 15 levels 
all of them being treated as single levels. The energy levels 
and oscillator strengths for the all transitions are taken from 
Wiese et al. (1966 1. The radiative probability of the hydro- 
gen 2s-ls two-photon transition and the distribution of the 
hydrogen two-photon continuum are given by Nussbaumer & 
Schmutz ( 1984). The photoionization cross sections of hydro- 
gen are computed with the formulae and tables of Karzas & 
Latter J1961> . The temperature-averaged free-free Gaunt fac- 
tor is interpolated from the table and accurate extrapolations 
given by Sutherland (.1998) . The electron collisional excitation 
and deexcitation rates for atomic hydrogen are calculated with 
the effective collision strengths taken from Scholz et al. ( I1990> . 
Callaway C1994) , Aggarwal et al. Cl991j , and Giovanardi et 
al. ( I1987> . The electron collisional ionization rates for atomic 
hydrogen are calculated using the approximate formulae of 
Johnson (J,972 1. The rate coefficients for charge transfer in col- 
Usions of with H and H^ with O are given by Stancil et al. 
lfT999ll . 

A.2. Singly ionized magnesium 

For singly ionized magnesium we use a model ion consisting 
of 9 single levels. The energy levels and oscillator strengths for 
the all transitions are compiled from the Atomic Spectroscopic 
Database (ASD) at the National Institute of Standards and 
Technology (NIST) and from Kurucz & Bell (IT995t . The pho- 
toionization cross sections of ionized magnesium are taken 
from the Opacity Project atomic database (Cunto et al. ll993> . 
The effective collision strengths for the electron excitation 
of ionized magnesium were computed by Sigut & Pradhan 
( 119951 . The electron collisional ionization rates for singly ion- 
ized magnesium are evaluated by fits of Voronov ( 1997), the 
approximate formula of Milkey & Mihalas (I1974> . and the 
semi-empirical formula of Lotz ( I1969> . 

A.3. Singly ionized iron 

For the specific case of the Fe II ion a 30-term, multilevel atom 
model is employed. A total number of 212 transitions between 
these terms from 2150 A to 55 jum are calculated. The energy 
levels and oscillator strengths for the all transitions in singly 
ionized iron are compiled from ASD at NIST and from Kurucz 
& Bell ( I1995> . The photoionization cross sections of ionized 
iron are taken from Nahar & Pradhan (1994j and Bautista & 
Pradhan ( I1998t . and averaged over resonance structures ac- 
cording to Bautista et al. ( I1998I I. The Maxwellian averaged 
rate coefficients for the infrared, optical and ultraviolet tran- 
sitions in ionized iron are taken from Zhang & Pradhan ( I1995> . 
The collision strengths for some transitions are calculated by 
using an empirical formula derived from the available experi- 
mental and theoretical data for Fe II (Li et al. ll993> . The rate 
of direct collisional ionization for singly ionized iron from the 
ground state is computed using the fits from Voronov (.19971 
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and the electron collisional ionization rates from excited levels 
are evaluated by the semi-empirical formula of Lotz (.1969j . 

A.4. Singly ionized barium 

For the case of the Ba II ion a 17-term, multilevel atom model 
is used with a total number of 56 transitions between these 
terms from 1500 A to 58 /im. The energy levels and oscilla- 
tor strengths for the all transitions in singly ionized barium are 
taken from Kurucz & Bell (.1995J . The photoionization cross 
sections of ionized barium are given by Butler ( 2000 1 and av- 
eraged over resonance structures according to Bautista et al. 
<I99 8 1. The electron collisional excitation rates for ionized bar- 
ium are taken from Schoning & Butler (1998 1 and compiled 
from Sobelman et al. ( 1981 1 and Van Regemorter ( 1962 1. The 
electron collisional ionization rates are evaluated by the semi- 
empirical formula of Lotz Cl969j . 

A. 5. Other atoms and ions 

The following elements are included in the non-steady study: 
He, C, N, O, Ne, Na, Mg, Si, S, Ar, Ca, Fe, and Ba. All elements 
are treated with the three ionization stages. The atoms and ions 
are assumed to consist of the ground state and continuum ex- 
cept the singly ionized elements with a detailed energy level 
structure. Atomic weights and ionization potentials are from 
Allen (Q973l. The partition functions are calculated with the 
polynomial approximation fit obtained by Irwin (19811. The 
photoionization cross sections of atoms and ions are evaluated 
with data of Verner & Yakovlev ( 1995 1 and Verner et al. ( 1996 1. 
As for H the temperature-averaged free-free Gaunt factor is 
taken from Sutherland ( 1998 1. The free-free absorption coef- 
ficient is calculated with the effective nuclear charge including 
screening effects (Sutherland & DoDita ,1993J . The electron col- 
lisional ionization rates for atoms and ions are computed using 
the approximate formulae of Voronov ( I1997> . 

Appendix B: Molecular data 

B. 1. Negative hydrogen ion 

The photoionization cross section data for negative hydrogen 
ion are taken from Wishart ( I1979> . The free-free absorption 
coefficient of negative hydrogen ion was calculated by Bell & 
Berrington ( fT987l . 

B.2. Other hydrogen molecules 

The neutral-neutral reactions for hydrogen molecules are listed 
in Table IB. II with their names and references from which 
the corresponding rate coefficients are taken. Tables IB. 21 and 
IB. 31 list the ion-neutral and ion-ion reactions for hydrogen 
molecules with the relevant references, respectively. The elec- 
tron reactions for hydrogen molecules are listed in Table IB.4I 
with their names and references from which the correspond- 
ing rate coefficients are taken. Table IB. 51 gives the photoion- 
ization and photodissociation processes, their names, and ref- 
erences from which the photoionization and photodissociation 



Table B.l. List of neutral-neutral reactions. 



iNO. 


Reaction 


Source 


1 


H + H-»H++H + e- 


Hollenbach & McKeeTW 


2 


H -1- H -» H+ + H- 


Hollenbach & McKee.1989. 


3 


H + H + H^H2+H 


Palla et al.[T983l 


4 


H + H + H2 ^ 2H2 


Palla et al.[T983l 


5 


H(2s) + H ^ H+ -H e- 


Rawlings et al.fTW^ 


6 


H(2p) -H H ^ H+ + e- 


Rawlings et al. 1993J 


7 


H(3) + H -> H+ + e- 


Rawlings et al. 199j] 


8 


H(2s) + H(ls) ^ H2 + r 


Latter & Black 19<^ 


9 


H(2p) + H(ls) ^ H2 + 7 


Latter & Blackll99n 


10 


H + H2 ^ 3H 


Roberge & DalgamofT9SII 


11 


H + H2 -» H + H+ 


Hollenbach & McKee 19W 


12 


H + H2 ^ H+ + H + e- 


Hollenbach & McKee.1989. 
Hollenbach & McKee[T989l 


13 


H + H2 ^ H+ + H2 + e- 


14 


H* + H2 ^ H+ + e- 


Culhane & McCray'T995] 


15 


H2 + H2 ^ H+ + H2 + e- 


Hollenbach & McKee.19891 


16 


H2 + H2 ^ H2 + 2H 


Rawlings 1988] 



Table B.2. List of ion-neutral reactions. 



No. 


Reaction 


Source 


17 


H- + H ^ 2H + e- 


Abel et al. 1997 


18 


H- + H ^ H2 + e- 


Galli & Palla.19981 


19 


H+ + H ^ H^ + 7 


Galli & PallafTWg] 


20 


H+ + H2 HJ + 7 


Galli & PallalWH] 


21 


H+ + H2 H+ + H 


Galli & Palla 1998] 


22 


H+ -1- H -» Hj H+ 


Galh & Palla 1998' 


23 


H+ + H2 ^ H+ + H 


Hollenbach & McKeeTW 


24 


H^ -1- H2 ^ H2 + H+ + H 


Hollenbach & McKee 1989] 


25 


H+ + H ^ H2 + H^ 


Galli & Palla 1998 


26 


HJ + H2 H2 + HJ + H 


Hollenbach & McKee[T989l 


27 


H+ + H2 ^ 2H2 + H+ 


Hollenbach & McKee[T989l 



Table B.3. List of ion-ion reactions. 



No. Reaction Source 

28 H- + H+ ^ 2H Galli & Palla 1998J 

29 H- + H+ H^ + e- Galli & Palla il998l 

30 H- + H+ ^ H + H2 Abel et al. fTWTI 



cross sections are taken. These cross sections are used in calcu- 
lating the ionization and dissociation rates in the radiation field 
presented in Sect. 13. II 

Appendix C: Rate equations for H+, H~, H2, HJ, 
and H+ 

Using the physical values and quantities from Sect. l3.5l the rate 
equation for ionized hydrogen may be written as 

dNu* ^ _Nw_ gY 

15 15 
k=\ k=2 
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Table B.4. List of electron reactions. 



No. 


Reaction 


Source 


31 


e + H- ^ H + 2e- 


Abel et al.TWT' 


32 


e" + Hz ^ 2H -1- e" 


Stibbe & TennvsonlT^^ 


33 


e- + Hz ^ H + H- 


Fuller & Couchmanl2000l 


34 


e- + Hz ^ H+ + 2e- 


Hollenbach & McKeefT989l 


35 


e- + HJ ^ 2H 


Galli & Palla 195^1 


36 


e- + H+ ^ H + Hz 


Galli & Palla 1998| 


37 


e- H+ ^ 3H 


Hollenbach & McKeefT989l 


Table B.5. List of photoionization and photodissociation pro- 


cesses. 




No. 


Process 


Source 


1 


Hz + 7 ^ H+ -H e- 


Yan et al.n^EDOTl 


2 


Hz + 7 ^ 2H 


Abel et al.lTWTl 


3 


Hz + 7 ^ H; ^ 2H 


Abel et al.[l997l 


4 


H+ + 7 -» H -H H+ 


Stancil f^^H 


5 


HJ + 7 2H+ + e- 


Shapiro & KangfTWTl 



15 15 
k=l k=l 

+ (ki + kz)Nl, + fcoA^HoA^H, 

+ kziNu^Nu" + k24NH*Nu2 + ^27A^h+A^H2 

+ koHNo*Niio + (R4 + 2R^)NHi 

- kigNH*Nu« - (kiQ + k2i)Nn*Nii, 

- (kis + k29)NH-NH* - knoNH^Noo , 



(C.l) 



where R4 and R5 are the photodissociation and dissociative 
photoionization rates evaluated with the relevant cross sections 
from Table IB. 51 respectively; the rate coefficients for the in- 
volved reactions are taken from Tables IB.1HB.3| Radiative re- 
combinations to the ground state of hydrogen lead to imme- 
diate reionizations from the ground state, so they are omitted. 
The rate equation for negative hydrogen ion is 



dt 



Nn- 



+ NiioNeaH- -Nh-Ph- 



+ k2Nlo + kuNnoNn. + ^bsA^cA^h, 

- (kn + kis)NH-Niio - kioNu-Nn* 

- (k2S + k29)Nn-Nn* - k^iNeNn- , 



(C.2) 



where au- and Pu- are the total coefficient of photo-attachment 
of neutral hydrogen and electron to negative hydrogen ion and 
the photo-detachment rate of negative hydrogen ion in the con- 
tinuum radiation field, respectively; the rate coefficients for the 
corresponding reactions are taken from Tables IB.1HB.4| The 
rate equation for molecular hydrogen H2 is given by 



dt 



+ kisNa-Niio + k22NH*Niio + fesA^H+A^H" 



+ k2iNH*Nii, + kioNn-Nn- + k^^N^Nii* 

- (kiQ+ku +ki2)Nu»Nn2 

- ki4NH2NH0 - (ki5 + he)N^,_ 

- (k2Q + k2l)Nii*NH2 - k2iNH*NH2 

- ihl + hi + ^34)A^eA^H2 , 



(C.3) 



where R\, R2, and Rj are the photoionization, direct photodis- 
sociation, and two-step photodissociation (the Solomon pro- 
cess) rates of molecular hydrogen H2 evaluated with the cross 
sections from Table IB.51 respectively; the rate coefficients for 
the involved reactions are taken from Tables IbTMB .41 Finally, 
the rate equations for molecular hydrogen and H3 are 



-^+R^Nn2-(R4+R5)Nn* 

(ksNiioas + hNu\2p + kiNH«,3)NH0 
(kn + ki2)NiioNii, + ^isA^Hj 

^19A^H-A^H» + ^2lA^H+A^H2 + k25Nii*NH» 
k22Nn*NiiO - (^23 + k24)NH*NH2 

kioNn-Nu* - kiiN^Nw 



(C.4) 



and 



■ kuNH.NH" 



exp 



(C.5) 



+ k2QNH*NH2 + k23NiiiNH2 

- k25NH*NiiO - (k26 + ^27)A^H(A^H2 

- (^36 + ^37)A?eA'H+ , 

respectively. 
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